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Abstract 
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general properties of their solutions and find the drift criterion for chain motion. For ratchet 
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I. INTRODUCTION 



During the last decade, much attention was devoted to the study of the so-called 
ratchet effect, i.e., the rectifying of nondirectional stochastic or deterministic driving forces 
into directional motion of (quasi-)particles, for comprehensive reviews, see, for example, 



Refs 



.mm. 



This phenomenon is o 



Drominent practical importance and constitutes 
3, particle segregation Q, reduction of a 
and many others 

m. 



IG 



a theoretical basis for Brownian motors 
™te. density H,—ngofac,vsta.„;e;.faceB 

There exists a large variety of systems that exhibit the ratchet effect. Those of them 
where transported particles are exposed to a driving force and move in an asymmetric pe- 
riodic (ratchet) potential form a class of so-called rocking ratchets In turn, stochastic 
and deterministic rocking ratchets are distinguished, depending on whether they are driven 
only deterministically, or whether they are subjected also to stochastic forces. The fiuctu- 

ib-threshold 



ating noise component in turn allows for activated escape events ^| even at su 
driving. The stochastic situation has been studied for both, noninteracting 12 



interacting 



3, 12, 18, Q 



14 



3 



and 



particles, and the strong influence of the inter-particle interac- 



tions was revealed. In contr ast, the deterministic case has been investigated primarily for 



noninteracting particles 



21 



22 



23 



2J]. At variance with earlier work [25| focusing 



on the transport of driven linear defects (i.e. elastic chains) which diffuse on asymmetric 
substrates at finite temperatures, the objective of this study is to investigate the transport 
properties of a chain of charged, interacting particles in deterministic rocking ratchets. This 
model can be used for the investigation of one- dimensional crystals in carbon nanotubes 3| 
and ion transport through synthetic nanopores 27 1. 

The paper is organized as follows. In Sec. II, we derive the reduced system of two 
equations that describes the deterministic, overdamped motion of a chain. We carry out the 
general analysis of this system in Sec. III. Specifically, we introduce a class of its steady-state 
solutions, analyze the chain dynamics in the inverted potential, and derive the drift criterion 
of a chain. In Sec. IV, we study analytically and numerically the overdamped transport of 
a chain in the double-sine ratchet potential. Concluding remarks are contained in Sec. V, 
and the phase-modulated sine potential is introduced in the Appendix. 
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II. EQUATIONS OF MOTION 



We study a chain of charged particles of identical masses M that interact among each 
other via the Coulomb interaction, a repulsive interaction, and additionally with a (sub- 
strate) ratchet potential V{x) having the period d. Moreover, the particles are driven by 
a longitudinal alternating electric field E{t) of a temporal period 2T. Therefore, the total 
potential energy of a chain V^ot includes the interaction energy \4nt and the potential energies 
Vr and Ve\ formed by a ratchet potential and an electric field, respectively. Assuming that 
any neighboring particles have opposite charges, i.e., gj+i = —qi and \qi\ = q, and a repulsive 
interaction depends on the inter-particle distance (t) is the coordinate 

of the ith particle, r > 1 since it prevents the unbounded contraction of a chain], these three 
parts of the total energy can be written as follows: 



O ' ^ />'■ ry> O ■ ^ ,~fi ,~fi IT 

^ jj I J-j J'j I ^ I J-j J-j I 

V, = Y.V{xi), V,, = -J2qiXiE{t). (2.1) 

i i 

Here 6 is a dimensioned parameter characterizing the strength of the repulsive interaction, 
and the primes on the summation signs imply that i ^ j- 

Along with the potential force —dVtot/dxi on the ith particle a friction force —Xxi (A is a 
damping coefficient) is also acting. Then the equation of motion of each, individual particle 
reads 

dV 

Mx, + Xx, + -^ = q,E(t) + /(x,), (2.2) 

where f{x) = —V'{x) (here and below, the prime denotes the derivative with respect to 
the argument of the function) is a force field generated by a ratchet potential. The set of 
equations ()2.2j) is very complicated for a detailed analysis of the chain motion because they 
are coupled and contain the nonlinear terms dVi^t/dxt and f{xi). Nevertheless, the problem 
of coupled ratchets does possess a very rich behavior, such as anomalous hysteresis, self- 



oscillations, absolute negative mobility, etc. 28|, and thus is demanding to analyze. Since 



our aim is to study a ratchet mechanism of the chain transport, the nonlinear nature of the 
ratchet force f{xi) must be taken into account, while the nonlinearity of the interaction force 
dVint/dxi seemingly is not so essential. Therefore, in order to partially simplify the problem, 
we restrict ourselves to a harmonic approximation for V^nt- This approximation is valid if 
during the period 2T of the action of the electric field E{t) the particle displacements are 



much less than the equihbrium distance a [at V{x) = and E{t) = 0] between the nearest 
particles. 



According to the findings in 29[, a periodic chain with equidistant particles exists only 
if r > r2 ~ 2.799; in the opposite case the minimum of Vi^t occurs for a chain with infi- 
nite period. Assuming that the condition r > r2 holds, it is convenient to introduce the 
coordinates 2na + and {2n + l)a + x~ [n {n = 0, ±1, ...) numbers the chain cells which 
contain two particles (see Fig. ^] of positive and negative charges, respectively. The equa- 
tions of motion for the displacements x^ = x^(t) and x^ = x^(t) of these particles from 
ihe equilibrium positions follow from Eq. ()2.2j) . Using the harmonic approximation for V^nt 
2^ and the notations E{t) = Eh{t) and f{x) = fQg{x), where E is the amplitude of E{t), 
/o = |min/(x)|, and h{t) and g{x) are the dimensionless electric and ratchet force fields, 
respectively, the equations of motion assume the form 

x+ + 2Qxx^ + H ^2(„_„)_i(x+ - X-) 

m 

+ B2(n-m){xl - x+) = Ah{t) + Rg{xl), 

X~ + 2Qxin +^^Y1 B2(n-m)+l (^^ - X'^) (2.3) 
m 

+ B2(n-m){Xn " X^;^) = -Ah{t) + Rg{x~). 

Here Vtx = A/2M, fi^ = q^Ml^ I = (6/g2)i/('-i) ig the length scale, A = qE/M, R = fo/M, 

(2.4) 



Br, 



-1)" ^ (r + l)ln2 



|nP C(^)l^l''^^ 



7i(r) = [r((r) / \n2Y^^^~^\ a = /7i(r), and C{r) = I]^!?^"'' is the Riemann zeta function. 
Note that Eqs. (2.3) are valid if the condition \x^ — x^\ <^ a holds for all n and m. It 
provides the applicability of the harmonic approximation and, as can be shown readily, it 
does not prohibit the existence of unbounded (at t — > oo) solutions of these equations that 
describe the chain transport. 

The system of equations (2.3) constitutes a theoretical basis for the study of ratchet- 
transport of ionic chains within the harmonic approximation for the inter-particle interac- 
tions. Unfortunately, at present there are no methods to find and analyze its unbounded 
solutions. Since this system contains an infinite number of coupled nonlinear equations, the 
development of such methods is a very complicated problem which, in general, can be solved 
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only approximately. However, we will show that, in a particular case, the chain transport 
can be studied in detail without approximations. 

More specifically, we consider the case when the equilibrium chain period a in absence 
of the ratchet potential V{x) is a whole integer of the potential period c/, i.e., a = Ld, 
where L is a natural number. Although this commensurability assumption is rather restric- 
tive (small discrepancy between a and Ld violates the chain periodicity in presence of the 
ratchet potential), it permits us to reduce the infinite system (2.3) to the system of only 
two equations. Indeed, if the initial conditions x^{0) = x± and x^(0) = v± {x± and v± are 
the initial displacements and velocities of the positive and negative charges) hold for all n, 
then all positively charged particles and all negatively charged particles move identically. 
Designating in this case = x^, from Eqs. (2.3) we get 

x+ + 2nxx+ + V(l)(x+ - X-) = Ah{t) + Rg{x^), 

(2.5) 

X- + 2nxx~ + icj2(l)(x- - x+) = -Ah{t) + Rg{x-), 

where t^^(l) = 2fi^ J2mB2m-i is the squared frequency of optical vibrations of chain parti- 
cles that corresponds to the dimensionless wave number k = 1 f2^]. Note that Eqs. ()2.5|1 
are derived from Eqs. (2.3) without approximations and they precisely describe the chain 
dynamics in the examined case. However, since a nonlinear oscillator driven by a periodic 
force can exhibit chaotic behavior that is characterized by strong sensitivity to initial con- 
ditions [3^, the chain dynamics governed by Eqs. (2.3) and ()2.5|) can be quite different if 
the corresponding initial conditions slightly differ. With increasing damping constant the 
chaotic domain in the parameter space is reduced, therefore we expect (and this is confirmed 
by simulations) that it vanishes in the overdamped limit {Qx 00). In other words, the 
chain dynamics in this limit is expected to be regular and predictable. 

Introducing the dimensionless time r = t/2T and the dimensionless particle displacements 

x^ {2Tt) / d, from Eqs. ()2.5|) we obtain in the 

\{w -u) + ^iG{w), 

(2.6) 

- \{u -w)+ fiG{u). 

Here x = ^\/Tuj'^{l) and = A/duj'^{l) are the dimensionless parameters characterizing the 
electric field frequency and amplitude, respectively, H{t) = /i(2Tr), /i = R/duj'^{l) is the 



w = w{t) = x'^{2Tt) / d and u = u{t) = 
overdamped limit 

= -mr) - 

dr 
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dimensionless parameter characterizing the ratchet force amphtude, G{x) = g{dx), x = x/d 
is the dimensionless coordinate, and according to Eq. ()2.4|) 



3 ^ 

r-l 



x^-7C(3, 



41n2(l-2-"-2)(r + l) 

(2.7) 



Assuming that all particles at t = are in equilibrium (their equilibrium positions are 
those in absence of the ratchet potential), the initial conditions for Eqs. ()2.fij) are written 
as w(0) = m(0) = 0. The system of equations (j2.6|) provides a very useful tool for studying 
the transport properties of ionic chains. Indeed, on the one hand, it is a rather simple set 
of only two coupled, ordinary differential equations of first order and, on the other hand, it 
accounts for the inter-particle interactions and the action of the ratchet force. 

By the definitions, the functions H{t) and G{x) have unit periods, H{t + 1) = H{t) 
and G{x + 1) = G(x), and zero mean values, /q^ H{T)dT = and Jq G{x)dx = 0. The latter 
condition shows that the total work, delivered by the ratchet force field f{x) on any interval 
of length d, equals zero. In general, the functions H{t) and G{x) can be both, continuous 
and discontinuous. But here, to simplify the numerical solution of Eqs. ()2.6|) . we consider 
them as smooth, differentiable functions. 

In what follows, to find the drift criterion of a chain, we shall use the dimensionless 
potential energy U = U{w,u,t) that reduces Eqs. ()2.(i|) to the form 

dw d ^ du d ^ 

dr dw ' dr du 

Introducing the representation V{x) = VoW{x), where Vq > and W{x) is the dimensionless 
ratchet potential, and by use of the definitions of f{x) and /o, we get /o = {Vo/d) ma.xW'{x) 
and 

G{x) = -W'{x)/ma.xW'{x). (2.9) 
Finally, from the inspection of the right-hand sides of Eqs. ()2.6p and ()2.8p we obtain 

U = — 0sin(27rr)(tt; — u) + ^{w — 

+l^[W{w) + W{u)]/maxW'{x). (2.10) 
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III. GENERAL RESULTS 
A. Periodicity analysis 

In order to exclude from consideration transient processes, we need to examine the asymp- 
totic, steady-state solutions of Eqs. ()2.6|) . These solutions depend on many factors, such as 
the form of a ratchet potential, characteristics of a chain, initial conditions, etc., and con- 
sequently can be studied in detail only numerically. However, using the periodicity of H{t) 
and G{x), it is possible to introduce different classes of the steady-state solutions. One par- 
ticular such class is generated by those periodic solutions of Eqs. (j2.6|) that asymptotically 
r ^ oo obey: 

w{t + k) = w{t) + K, u{t + k) = u{t) + K, (3.1) 

where k and K are natural and integer numbers, respectively, that have no common factors. 
In this case, the periodicity and increment of the functions w{t) and u{t) are described by 
the pair {k, K} which at K ^ corresponds to the drift state of a chain. Since the reduced 
chain displacement Aw = \imr^oo[w{T + k) — w(r)]/A; that occurs during one period of H{t) 
is given by Aw = K/k (note that the harmonic approximation is valid if \K\/k <^ L), we 
shall term the chain transport characterized by the pair {k,K} as "integer" if A; = 1, and 
"fractional" if > 2. 

According to the conditions (jH.lj) . the average velocity of a chain or drift velocity v = 

((i/2T) limT-^oo w;(r)/r [= ((i/2T) hmT-^oo w(r)/r] is reduced to 

d ,. w(t + k) — w(t) ,„ ^, 

t, = _ lim ^ 1 ^ 3.2 

which, in turn, yields v = VqV, where Vq = duj'^{l)/2Qx and v = xK/k is the dimensionless 
drift velocity. Taking into account that the periodicity and drift parameters k and K depend, 
in general, on all parameters of Eqs. ()2.6|) . we conclude that w is a discontinuous linear 
function of x- 

We emphasize that this class does not exhaust all the steady-state solutions of Eqs. ()2.6p . 
Moreover, the symmetry approach does not permit us to find their unique steady-state 
solution in each concrete case. Therefore, to study the transport properties of a chain, it is 
necessary to numerically find the solution of Eqs. ()2.(jj) (with zero-valued initial conditions) 
and examine its long-time behavior depending on the form of a ratchet potential, the electric 
field characteristics, and the chain parameters. 
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B. Chain dynamics in the inverted potential 



We now consider the chain dynamics in the inverted potential Vin{x) = V{—x) that 
generates the reduced force field G'in(a;). According to Eqs. (j2.6|) . in such a potential the 
displacements Wi^^r) and Ui^^r) of positively and negatively charged particles from their 
equilibrium positions are governed by the equations of motion 

^dWin ^ (j)H{T) - ^{Win - Win) + /iGin(Win), 
X-J^ = -(t>H{T) - i(Min - Win) + /iGin(Min) 

[win(O) = Min(O) = 0]. Taking into account that G'in(— x) = — G'(x), one can show from 
Eqs. fl2.6p and ()3.3|) that Wi^ij) = —u{t) and Ui^ij) = —w{t). This implies that if the chain 
dynamics is known for the ratchet potential V{x), then it is also known for the corresponding 
inverted potential Via{x) as well. In particular, if a chain in the potential V{x) drifts along 
the X-axis, then in the inverted potential Vin{x) it drifts in the opposite direction with the 
same average velocity, i.e., fin = —v. We emphasize that in refiection-symmetric potentials 
the drift state of a chain, i.e., a chain state with v ^ 0, does not exist. Indeed, as it is shown 
above, the general condition = —v must hold. On the other hand, if V{—x) = V{x) 
then Vij^{x) = V{x) and so the condition t>in = v also must hold. It is obvious that both 
conditions are met simultaneously only if f = 0. 

The above mentioned features of the chain djTiamics permit us to study the overdamped 
transport of ionic chains only in those ratchet potentials that induce the chain drift, say, with 
positive velocity v. In the following we consider the simple-structured ratchet potentials. 
We assume that the reduced potentials W{x) and the corresponding force fields G{x) have 
only one maximum and one minimum on unit period and, in addition, maxG'(x) > 1, see 
Fig.121 

C. Drift criterion 

To find the conditions that lead to the drift state of a chain, we rewrite Eqs. (I2.fi|) as 

+ u) = iJ,[G{w) + G{u)], 
8 



where the second equation is obtained by summing ()2.6|) . If /i/x ^ 0, then it reduces to 
the equation d{w + u)/dT = 0, which with w{0) = m(0) = yields w{t) = —u{t). Using 
this relation and the condition that /x/x — > 0, the first equation in (|3.4|) takes the form 
xdw/dr + w = (pH^r). Its solution 



^ r H{T-T')e-^'/^dT' (3.5) 



W[T) = 

X 

shows that 

w(t + 1) = w(t) + ^e""/^ [' if(-r')e'"'/^c/r', (3.6) 
X JO 

and so w{t+ 1) = w{t) for r oo. This means that the drift state of a chain does not exist 
if fi/x 0. Since the parameter x is proportional to the electric field frequency, its decrease 
leads to an increase of the maximal particle displacements, yielding f |;^>o = if v\^=o = 0. 
The last condition is violated if the amplitude parameter is large enough. Therefore, to 
find the drift criterion of a chain, we need to consider its dynamics as x 0. 

According to Eqs. ()3.4|) . in the stationary regime {x 0) the chain dynamics is described 
by the system of nonlinear equations 

(PH{t) - -u)+ iiG{w) = 0, 

(3.7) 

G{w)+Giu) =0. 

If for each time the chain energy U has a minimum value, i.e., if in virtue of Eqs. ()3.7|1 the 
conditions 

^>0 ^^-(^R-Wo (38) 
dw"^ ' dw'^ du^ \dwdu J 

hold, then v\^=o = 0. The latter condition in (jSD), G'{w) + G'{u) - 2fiG'{w)G'{u) < 0, is 
weaker than the former, 2fiG'{w) — 1 < 0. Hence, it is violated first under increasing of the 
parameter 0. Let 0cr be the critical value of the parameter (p such that, for (p = (per, the 
latter inequality in ()3.8|) is reduced to equality at some instant of time. It is obvious that 
this occurs for the first time at r = ri, where ri(< 1) is the minimal solution of the equation 
H{ti) = 1. Then, using the first equation in ()3.7|) . we obtain 

0cr = liw^r - u~) - f^G{w-), (3.9) 
where w~ = w{ti — 0) and u^^ = u{ti — 0) are defined by the system of equations 

G'{w-,) + G'(n-J - 2^^G'{w-^G'{u-;) = 0, 

(3.10) 

G{w~) + G{u-,) = 0. 
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At = 0cr and r = ri + the chain particles instantly move to the new equilibrium 
positions = w{ti + 0) and u'^^ = u{ti + 0), which are defined by another system of 
equations 

0cr - 1{W+ - U+) + fiG{w+) = 0, 

(3.11) 

+ = 0. 

As Fig. 121 illustrates, in this case the positively charged particles pass into the next potential 
wells, while the negatively charged particles do not leave their own wells. A detailed analysis 
shows that during the second and each following periods of H{t) both types of particles 
instantly move into the next potential wells twice. In other words, during each period of an 
alternating electric field a chain in the steady-state regime is displaced by two periods of a 
ratchet potential. 

Thus, the drift criterion of a chain, that leads to the condition f |;^=o 7^ 0, has the form 
> 0cr- As — s> 0cr and X ^ a chain exhibits the integer transport with the drift 
parameter K = 2 and with an average velocity v = 2%. According to Eqs. ()3.9|1 and (jH.lOj] . 
to calculate 0cr it is necessary to know the explicit form of a ratchet potential. However, 
taking into account that for w = Xi and u = — 1 the conditions G{w) + G{u) = and 
v\^=o = hold, we find the general condition 

0cr > - X5 + 1) + /i, (3.12) 

which can be used for approximate estimation of 0cr- 

Note also that, because for noninteracting particles the chain energy ()2.10|) does not 
contain the term {w — m)^/4, the drift criterion of free particles assumes the form > /i. 
The main distinction between the drift states of interacting and noninteracting particles 
thus is the result that in the latter case — 5> oo as x — > 0. 



IV. DOUBLE-SINE POTENTIAL 



A. Analytical results 



We examine the chain dynamics in the asymmetric ratchet potential composed of two 



spatial harmonics 



12j |. This so-called double-sine potential is defined as V^^x) = VodWd{x), 
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where Vod is a positive constant and 



Wd{x) = - sin[27r(x + Xdj] - r/d sin[47r(x + Xd)]. (4.1) 

Here, the parameter r7d(> 0) characterizes the form of the reduced potential Wd{x), and 
the parameter Xd = Xd/d defines the positions of its extrema. The function Wd{x) and 
the corresponding reduced force field Gd{x) = —W^{x)/maxW^{x) [fd{x) = fodGd{x), 
fod = (Vod/d) ma.xW^{x)] both have qualitatively the same forms as those depicted in Fig. |2l 
i.e., they have only two extrema per unit period and mmWd{x) = Wd{0), if rjd < //od = 1/8 
and 



- 1 \/l + 32^d-l .... 

Xd = — arccos -^^ . (4.2) 

27r 87]d ^ ^ 

(Yet another ratchet potential that possesses the same properties as Vd(x) is introduced in 
the Appendix.) For these conditions, we find that ma.xW^{x) = 27r(l — 2?7d), m.inW^{x) = 
-27r(l + 2r]d), fod = 27r(l - 2r/d)K)d/c?, 

Gd{x) = {cos[27r(x + Xd)] + 2?7d cos[47r(x + Xd)]} 

x{l-27]d)-\ (4.3) 

G'd(xi) = minGd(x) = —1 if xi = 1/2 — Xd, G'd(x2) = G'd(O) = if X2 = 1 — 2xd, 
Gd{xi) = maxG'd(a^) = (1 + 277d)/(l — 2?7d) if X4 = 1 — Xd, and Gdix^) = Gd{xrj) = 1 if 

. - 1 VI + ISr^d - 1 
X3,5 = 1 - Xd =F 7^ arccos , (4.4) 

where the upper and lower signs correspond to the indexes 3 and 5, respectively. 

To calculate 0cr, we proceed as follows. First, instead of the second equation in ()3.1()j) . 
we introduce the two equations G{w^j.) = —p and G{u^j.) = p (—1 < p < 1). Then, taking 
into account the conditions xi < w~ < X3 and X5 — 1 < -u^jj, < xi, we find their solutions 

i±i-.,l..„os£(M^. (4.5) 

WrJ 2 27r 8r]d 



where Z{p) = y 1 — 16r]dp + 32?7j(l + p), and, using Eqs. ()4.3|1 and ()4.5|1 . we reduce the first 
equation in ()3.10|) to the form 

E r^'^^' ~ '^'^^ = ^p (4.6) 

r Z{ap)J6Ar^l-[Z{ap)~l]^ 2 
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((J = ±1). Since the left-hand side of Eq. ()4.6|) . L{p), is a monotonic odd function and 
L{p) —>■ oo as p ^ 1, this equation always has a unique solution with respect to p. If that 
solution is known, then from Eqs. (|3.9|) and (j4.5p we get the desired formula 

</'cr = 2 + /^P - ^ arccos — . (4.7) 

According to Eqs. ()4.(i|l and ()4.7j) . is a universal function of p and 1]^. A corresponding 
3-D plot, obtained via the numerical solution of Eq. (|4.6|) . shows, see Fig. |21 that 0cr is an 
almost linear function of these variables. 

An analytical solution of Eq. ()4.6p is possible only in some special cases. Specifically, if 
T^d — > then, calculating the leading term of L{p), Eq. ()4.6|) gives 

P(3 - 2p') p 
(1 - p^y/^ r]d 

At /i/r/d oo its approximate solution reads p = 1 — (rjd/irp)'^^^ /2, and Eq. (j4.7j) yields 
(per = 1/4: + p. Another example corresponds to the limit p —>■ 0. Because L{0) = 0, the 
solution of Eq. ()4.6|) tends to zero as /i — > 0, and so 0cr = 1/2 ~ Xd- We emphasize that the 
cases r/d = and p = are degenerate. This means that for r/d = and p = the drift state 
of a chain never exists, while for r]d ^ and p 7^ it is realizable. The reason lies in the 
breakdown of the spatial symmetry. 



B. Numerical results 



We solved Eqs. (12. 6p with zero initial conditions by the fourth-order Runge-Kutta method 
for H{t) = sin(27rr) and G{x) = Gd{x). The analysis shows that each steady-state solution 
of Eqs. ()2.6p satisfies the conditions ()3.1|) . At > the typical dependence of the 
dimensionless chain displacement Aw on x is depicted in Fig. |3] a. The changes of Aw 
occur in a very narrow intervals {xhxi) (^^e index s labels these intervals) of the x-axis. 
We found that if the parameter x does not belong to these intervals, then k = 1, and the 
chain dynamics is characterized by the pairs {1, K}. According to our terminology, a chain 
exhibits an integer transport if 7^ 7^ 0. Its main features are as follows. First, the chain 
velocity v = Kx is a piecewise linear function of x that has a number of local maxima, 
see Fig. 0] b. These maxima occur due to the discrete character of the drift parameter K. 
Second, K\^=q is an increasing, step-like function of (p that equals zero if < (per and takes 
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TABLE I: The numbers k and K vs x for x ^ iXi^xi)- 



A/ 

A 


h 




X 


k K 




LOO 




n 479 


9 1 


0.4693878 


39 


38 


0.473 


9 2 


0.469388 


21 


20 


0.4732 


8 1 


0.46939 


7 


6 


0.47327 


47 2 


0.4694 


3 


2 


0.473278 


82 1 


0.47 


3 


2 


0.4732785 


150 1 



even numbers if > 0cr- Specifically, in accordance with the analytical results, K\^=q = 2 at 
~ 0cr- To illustrate for x << 1 the chain dynamics in the drift state, the time dependence 
of the particle displacements w{t) and u{t) and the phase trajectory of the chain motion 
are shown in Figs. El and El for = 3 and fi = 2. Since a chain reaches the steady state only 
if r >> X (the relaxation time is of the order of x); the functions w{t) and u{t) represent 
the steady-state solution of Eqs. ()2.(i|l . which satisfy the conditions w{r + 1) = w{r) +4 and 
u{t + 1) = u{t) + 4, at r > 1. Finally, the increase of x leads to the stepwise decrease of K 
and a smoothing of w{t) and u{t). 

Solving Eqs. ()2.6|) for x ^ {Xii X2)' discovered that the chain dynamics is characterized 
by the pairs {/c, K} with k > 2 and K ^ 0. In this case a fractional transport of a chain 
is realized. The number of the intervals (Xi;X2) equals K\^=q (s = 1, 2, i^|^=o) and 
their width grows with s. Within each such interval the chain displacement Aw = K/k 
assumes a stepwise function of x that satisfies the condition — 1 < K/k < Kg, where 
Kg = K\^=Q — s + 1 is the drift parameter to the left of the interval {xl, X2)- X approaches 
its boundaries on the inside, then k is strongly increased, K/k ^ Kg as x ^ Xh and 
K/k ^ Kg — 1 a.s X ~^ Xi- Table 1 illustrates these properties for the case represented in 
Fig. Hat X e ixtxt), where xl ~ 0.4693877 and xt ~ 0.4732787. To illustrate the chain 
dynamics at x G {Xi^X2)^ the time-dependence of the particle displacements w{t) and u{t) 
and the phase trajectory of the chain motion are depicted in Figs, d and |Hl respectively. If 
X > X2 then a drift of the chain does not exist and the phase trajectory explores a finite 
region of the phase plane. 

The fractional transport occurs also for noninteracting particles. But the x-intervals. 



13 



where such a transport exists, are much more narrow as compared to those observed in the 
interacting case. In particular, according to the above results X2 ~ Xi ~ 3.9 x 10^'^, while 
for noninteracting particles the width of the corresponding interval approximately equals 
1.4 X 10"! 

Note that in the phase-modulated sine potential, introduced in the Appendix, the ionic 
chain exhibits qualitatively the same features. We expect, therefore, that the results of the 
overdamped ionic chain dynamics are typical and robust for the considered class of ratchet 
potentials. 

V. CONCLUSIONS 

We have investigated the overdamped transport of a chain of charged, interacting par- 
ticles driven by a longitudinal alternating electric field that additionally interact with a 
smooth, non-symmetric, periodic ratchet potential. Assuming that the equilibrium parti- 
cle positions coincide with the minima of a ratchet potential, we have reduced the infinite 
system of equations that describes the dynamics of each chain to the system of two equa- 
tions which effectively describe the dynamics of only two, positively and negatively charged, 
representative particles. The reduced system of equations ()2.6p has the advantage of being 
particularly simple because it consists of two ordinary differential equations of first order, 
which are driven by the external force. 

Using the time-periodicity of an alternating electric field and the space-periodicity of a 
ratchet potential, we have introduced a wide class of corresponding steady-state solutions 
of Eqs. ()2.6|) . The mathematical structure corresponds to the drift state of the ionic chain. 
Particularly, all the steady-state solutions of Eqs. ()2.6|) obtained numerically in the cases 
of the double-sine and phase-modulated sine potentials belong to this class. Studying the 
chain dynamics in the original ratchet potential and in its inverted realization, we have 
shown that, depending on the parameter regime, a chain either does not drift at all for both 
realizations, or it has a finite drift velocity v which is opposite in value, and correspondingly 
— f , for the inverted ratchet potential. Considering the chain dynamics in the stationary 
regime, we have derived the drift criterion of a chain. Accordingly, the drift state of a chain 
takes place if during the first half-period of an alternating electric field the chain particles 
perform stick-slip transitions. 
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Applying analytical and computational methods for analysis of the chain dynamics in 
the double-sine and phase-modulated sine ratchet potentials, we have shown that the chain 
displacement, which occurs during one full period of an alternating electric field, is a mono- 
tonically decreasing, stepwise function of the electric field frequency. This function, scaled 
by the ratchet potential period, takes on only integer and fractional values. Therefore, only 
two types of the chain transport, namely integer and fractional, do exist. Both types oc- 
cur for tailored frequency intervals; the frequency intervals, however, that correspond to 
the fractional transport are a much more narrow than those corresponding to the integer 
transport. 
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APPENDIX: PHASE-MODULATED SINE POTENTIAL 

We define a new ratchet potential Vp{x), which we call the phase-modulated sine potential, 
as Vp{x) = VQpWp{x), where Vop is a positive constant and 



The dimensionless potential Wp{x) and the corresponding dimensionless force Gp{x) — 
—Wp{x) / maxWp{x) [fp{x) = fopGp{x), /op = (Vop/d) maxWp{x)] both have only two ex- 
trema on the unit period if the phase amplitude ?7p(> 0) satisfies the condition rjp < rj^p ^ 
0.31767, and minH^p(x) = M^p(O) if the parameter Xp [xp e (0,1/4)] is a solution of the 
equation 



In this case, ma.xWp{x) = 27r(l — rjp), mmWp{x) = — 27r(l -|- rjp), /op = 27r(l — r)p)VQp/d, 



Wp{x) 



sin{27r(x + Xp) + ?7p sin[27r(x + Xp)]}. 



(A.l) 



27rXp -|- r7pSin(27rXp) = 7r/2. 



(A.2) 




cos {2n{x + Xp) + rjp sin[27r(a; + a;p)]} 



x{l + rjpC0s[2TT{x + Xp)]}{l - rjp) 



-1 



(A.3) 
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Gp{xi) — —1 if xi — 1/2 — Xp, Gp{x2) — Gp{0) = if X2 = 1 — 2xp, and Gp{x4) — 
maxG'p(x) = (1 + ?7p)/(l - Vp) if = 1 - Xp. 
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Denisov, Denisova, Hanggi: Fig. 1 
FIG. 1: Schematic representation of the ionic chain. 
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Denisov, Denisova, Hanggi: Fig. 2 

FIG. 2: Form of the reduced ratchet potential W{x\ part (a), and the corresponding force field 
part (b), that are under consideration in this work. 
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Denisov, Denisova, Hanggi: Fig. 3 
FIG. 3: 3-D plot of (per as a function of fi and r/d- 



20 




Denisov, Denisova, Hanggi: Fig. 4 

FIG. 4: Plots of Aw, part (a), and v, part (b), vs x for = 3 and = 2. The intervals (xiiXl) 
(s = 1, ...,4) are not visible on this scale. 
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Denisov, Denisova, Hanggi: Fig. 5 
FIG. 5: Plots of the particle displacements w{t) and u{t) iov 4> = 3, h = 2, and x = 10 
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Denisov, Denisova, Hanggi: Fig. 6 
FIG. 6: Phase trajectory of the chain motion for the same parameters as in Fig. |21 
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Denisov, Denisova, Hanggi: Fig. 7 
FIG. 7: Plots of the particle displacements w{t) and n(r) for </> = 3, = 2, and x = 0.473. 
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Denisov, Denisova, Hanggi: Fig. 8 
FIG. 8: Phase trajectory of the chain motion for the same parameters as in Fig. |7| 
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